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Abstract 

Since the discovery of quasars in papers often appeared and appear the assertions that the redshift 
quasar distribution includes a periodic component with the period Az = 0.063 or 0.11. A statement of 
such kind, if it is correct, may manifest the existence of a far order in quasar distribution in cosmological 
time, that might lead to a fundamental revision all the cosmological paradigm. In the present time there 
is a unique opportunity to check this statement with a high precision, using the rich statictics of 2dF and 
SDSS catalogues (~ 85000 quasars). Our analysis indicates that the periodic component in distribution 
of quasar redshifts is absent at high confidence level. 

Keywords: (cosmology:) large-scale structure of Universe, (galaxies:) quasars: general, catalogues, 
methods: data analysis. 



1 Introduction 

As early as the first hundred galaxies with active nuclei and quasars have been discovered, the attempts to 
reveal the periodicity in their redshift distribution have been made. For example, the presence of peaks at 
z* = 0.061 • ?i, where n is the integer, for the dist ribution of 73 objects with non-thermal optica l continuum 
and z < 0.6 was mentioned in Burbidges papers (jBurbidge fc Burbidgd . Il957t iBurbidgd . Il968l ). They used 
this fact to confirm their hypothesis concerning the non-cosmological origin of the lines redshift in active 
galaxy and quasar spectra. However, the other interpretations are also possible. Particularly, some authors 
discussed the effect of the influence of occurence of several strong emission lines typical for quasars (Mg II, 
2800 A; C III, 1900 A; C IV, 1550 A; Ly Q , 1216 A), in the range of spe ctral observations (A > 3 300 A ), 
which might emulate the "humps" in N q {z) distributio n (see, for examp l e iKaritskava fc Komberd (1970l)). 

In t he fo l lowing years in a numbre of p apers (e.g. IJaakkolal d l 97lf): iTifftl (Il976l . Il989h: i Tifft fc Cockel 
(|l993f) : lTifftl (|l996l ); iBell fc Comeanl (f2TO3T ): lNarlikar fe Arpl (|l993h : iBell et "aTT (|2004 l) the authors reported 
about the observed quantization of the redshifts in spectra of near S-galaxies, which satisfies the "Tifft 
series" : 

Av r = 2-(°+ B / 9 \ (1) 

where v r is the radial velocity ex pressed in light velocity units, D and B are the integers. 

In iKhodvachikh (|l979Lll990h papers the cyclical changes of statisticslly brightest quasars were revealed 
in V filter, using the argument 

a; = In (1 + 2) (2) 

with Aa; « 0.19 period. In IKhodvachikh! (|1988l ) paper the dependence of a numbre of powerful radiopulsars 
has been plotted against x variable. In the centimeter wavelength region the cyclic changes were revealed 
with periods 0.12, 0.19 and 0.38. It is interesting that much later in iRvabnikov et al.l (|2001al lbl) papers, the 
similar result was obtained after the analysis of the redshist distribution of ~ 800 absorption lines in spectra 
of bright AGN with 0.10 < z < 3.7. The authors mentioned that a separate analysis of N a t, s (z) distribution 
in different celestial hemispheres indicates that the phase of periodicity is conserved. They interpreted such 
unexpected result as a sign of existence of an oscillatory regime in the Universe expansion and then as a 
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pr esence of a large scale ce l lular structure in the distribution of the absorption systems in quasar spectra. 
In lRvabnikov fe Kaminken (|2010l ) paper the presence of the periodic structure is considered on the basis of 

the a bsorbtion l ines. 

In lKarlssonl (|l97llfl97l papers the pe aks in N n (z) distribu t ion w ith t he step Ax ~ 0.19 w e re con sidered 
on the basis of 574 quasar data. In papers iBurbidge fc Napier ( 2001 ) and Bell fc McDiarmidl ( 20061 ) on the 
basis of a gross sample the peaks in N q (z) distribution were ment i oned at z* = 0.062, 0.3, 0.6, 0.96, 1.41, 



1.96, 2.63, 3.45. These peak values satisfied the "Bell series", [Belj (|200j): 



= 0.062 • (lOiV - M) 



(3) 



where N and M take on integer values (see the corresponding table in lBelll ( 20021 )). The initial series redshift 
coincides with the first term of the main Tifft series with D — 4, B — 0, that corresponds to the velocity 
« 18600 km/s. 

It is obvious that in order to confirm such a non-standard conclusion concerning the distinctive features in 
redshift quasar distrib ution it is necessary t o analyse a much wider statistica l information, which is included 



in the catalogues 2dF ( Croom et al. . 20041 ) and SDSS ( Shnider et all 2003). An d the papers with analysis 



2005) paper there was reported 



of such kind have really appeared. Thus, for example, in iTang fc Zhang 
the result of analysis of 290 quasars sample (the same sample as used in Karlsson and Burbidge papers) and 
the periodicity with the step Ax — 0.081 has been confirmed at 3 a level. However , the analysis o f large r 
samples (22497 quasars in 2dF and 46420 in SDSS-5), presented in the same paper ITang fc Zhangl (|2005l) . 
did not reveal any periodicity. From this they drew a conclusion that the point is in non-homogeneity of 
small samples and the selection effects, which the small samples are exposed. So, the problem, seemingly, 
can be abandoned, 



However, there is a numbre of papers where the authors, analyzing the large catalogues, find th e argu 



ments for the existence of the periodicity in N q (z) distribution. For example, in iBell fc McDiarmidl ([2006) 
paper, according to SDSS data 6 peaks in the power spectrum were d e tected with the st ep Az = 0.65 at z* = 
1.2, 1.8, 2.4, 3.1 and 3.7. And in papers lHartnett fc Hiranol (|2007al ) . iHartnettJ (|2007bl ) on the basis of 2dF 
and SDSS catalogues the authors reported the presence of the "humps" through Az = 0.0102, 0.0246, 0.0448, 
that assuming that Hq = 72 km/c Mpc corresponds to the cells of 44, 102 and 176 Mpc. Except that, th e 



distribution N a ( z) is represented in th e form z» = 0.062 • n, where n = 3, 4, 5, 6, 10, 20 (jHartnetti . l2007d ). 
Nevertheless, in ITang fc Zhangl (120051) paper the arguments have been expressed for the fact that SDSS 
catalogue may include a periodic component because of the selection effects. 

However, it is clear that the analysis of different quasar samples leads the authors to unlike conclusions, 
concerning the existence or non-existence of a far order in N q {x) distribution. At present there is a unique 
possibility to check with high precision the hypothesis about the possible periodicity, using the rich statistics 
of 2dF and SDSS catalogues (~ 85 thousand quasars). So, it is worth to consider this problem more 
accurately, and this is the goal of the paper. 



2 Observational data and periodicity extraction methods 



In th e current paper we used two catalogues: 2dF (jCroom et al.l . l2004l ) (22272 objects) and SDSS (jShnider et al 
20031) (63255 objects, release 7). Covering the celestial sphere by these catalogues is shown in Fig. Q] As 
one can see, these regions are overlapped and the density of the objects in 2dF catalogue is as much as one 
order larger than in SDSS. In last (seventh) release of SDSS catalogue some "gaps" have been filled with 
respect to the previous one, and as a result of it the sample become more homogeneous. 

The periodicity criterion is a function K(T) of the data and a trial period T, which takes on "large" 
values when the data includes the periodic component with period T, otherwise its values are "small". 

To investigate the periodicity we used four different criteria. The first of them and, apparently, the most 
familiar one is the Rayleigh criterion: 
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Figure 1: Covering the celestial sphere by 2dF (dark grey) and SDSS (light grey) catalogues. 

where xi are the values of the variable x from ^ for each quasar with redshift z, N - total amount of 
objects, and T - trial period in the units of x. Three other criteria analyse the structure of assumed periodic 
component (for variable stars it is called a light curve). For that purpose we subdivide the assumed trial 
period T in m parts (each of T/m length), calculate the phase of each quasar and add the unity in the 
appropriate m-th part of the period. As a result we obtain a histogram, which indicates the numbre of 
quasars that drop in each of m parts of assumed period T. Three criteria for that histogram analysis are 
the variety of the epoch superposition method and can be written as: 
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(7) 



where m is the number of parts (bins, intervals) in which the trial period is subdivided and rij (j 

1.2 in - the numbre of quasars which drop in the appropriate j-th part of a trial period. If the 

periodic component in the data is absent, then all m bins should have approximately the same values. On 
the contrary, if the data contain a periodic component they should strongly differ. All criteria have different 
sensitivity and reveal different aspects of periodicity, therefore it is not unreasonable to use all of them for 
a more reliable detection of the periodicity. 

If we use the stochastic data, the mean criteria values are: MK\ — 2, MK^ = 9, MK4 = 20, that 
corresponds to the accepted value m = 10, and the dispersions are, respectively, DK\ = 4, DK2 = 18, 
DK4 = 120. It means that if the criterion K\ takes on the value K\ = 7, then it exeeds the random signal 
level by 2.5 a. To calculate the theoretical values of MK3 and DK3 we s hould make u s e the formulas from 
Gurin et al. (1988). Statistical properties of the criteria are considered in iGurin et al.l (|1992h paper. 
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Figure 2: Quasar distribution over redshift in the SDSS catalogue. The absciss axis corresponds to the 
quantity x = ln(l + z) from @. The ordinate axis shows the numbre of objects in the x-direction interval 
of 0.004 width. 

3 Analysis of quasar distribution over redshift 

The quasar distribution over redshift z in two catalogues shown in Fig. [2] and |3] Analysing the plots one 
can draw a conclusion that the SDSS catalogue is more representative in the region of large and small x 
(x < 0.4 and x > 1.3). Four large maxima in the SDSS quasar distribution near w 0.5,0.8, 1 1.4 can be 
explained by the selection effect (the redshift z is detected using four spectral lines) and are not related to 
a periodicity. In general, we may observe in the plots some oscillations with smaller periods, which might 
appear a weak periodic component after a detailed analysis. 

The result of application of criteria (J4j) — {(7l) to the SDSS catalog in the interval < T < 0.14 is presented 
in Fig. 01 As it follows from the plots, all the criteria yield the similar results with diffrent confidence degree, 
therefore we use below the Rayleigh criterion only. The spectrum of the 2dF catalogue shown in Fig. [5] 

Indeed, in the SDSS spectrum there are maxima, mentioned by some authors, near the values T = Aa; = 

0. 063 and 0.11, though at a low confidence level. Except that there is one more, even higher maximum at 
T = 0.035. In the 2dF spectrum one can see the maxima at T w 0.034, 0.083 and 0.11, but the maximum 
near 0.064 is very low. Using approach of such a kind we cannot define exact values of the confidence levels, 
because for this purpose we have to simulate a stochastic sample with the same parameters as the one in 
Fig. [2] This approach, however, allows us to detect position of the periodic component with highest possible 
precision. Note that the values of the periodicity criterion in this case are very small. It is clear enough, 
because we try to extract a weak periodic component against a background of a very strong continuous 
signal, similar to investigating of behaviour of sea waves when measuring depth of an ocean. 

To increase the extraction reliability one should use another way. Namely: one should cut the background 
component and consider only the oscillations of AJV,(a;) with respect to a background. One can do it, 
subdividing the sample in narrow z-intervals and analysing the numbre of quasars in these narrow intervals, 

1. e. essentially, averaging the distribution inside each interval. The distributions in Fig. [2] and [3] are 
prepared using this technique. As a "background" we consider the mean value of 5 neighbouring intervals, 
where the interval of interest is in the middle. In all cases we consider the absciss as the middle of the 
appropriate interval. The result of application of this procedure to the SDSS catalogue are presented in 
Fig. [6] (for the 2dF catalogue the result looks the same) . If the periodic component in the distribution does 
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Figure 3: Quasar distribution over redshift in the 2dF catalogue. The ordinate axis shows the numbre of 
objects in the x-direction interval of 0.0046 width. 
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Figure 4: Spectrum K(T) from ((4]) - ([7]) of quasar distribution N q (T) according to the SDSS catalogue. 
The absciss corresponds to the trial period T . The criterion numbre marked in each panel. 
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Figure 5: Spectrum K\(T) of the quasar distribution according to the 2dF catalogue. Absciss corrsponds 
to a corresponding period T over coordinate x from (j2J . Rayleigh criterion (j4]) used. 

exist it should be distinguished in the plot even by naked eye. However the plot in Fig. |6] hardly looks like 
a periodic one and rather looks like a noise. It confirm s the spectrum of t he plot, calculated according to 



Rayleigh crinterion and shown in Fig.[7J As it is known (jGurin et al.l . ll992f) the mean value of the Rayleigh 
criterion for the stochastic data equals to 2 and the mean standard deviation is also 2. Indeed, we can 
reveal in Fig. [7] two maxima at = 0.063 and 0.111, however the confidence level does not seem to be 
high. One can only mark that these components are slightly stand out against their neighbours, but the 
reliable detection of the periodicity cannot be confirmed. Thus, according to the available data we cannot 
draw a conclusion that the periodic component in ln(l + z) coordinate presents in the quasar distribution 
over redshift. 

It, however, does not mean that the periodicity is absent for other coordinate choice as well. It is possible 
to check the periodic properties for other variables. The most interesting variable here is 



R(z) 



dz 



^n m (i + z) 3 + n A 



(8) 



which has the physical meaning of the geodesic cosmological distance. Note that when z — > oo R(z) tends 
to a constant value, i.e. to the horizon. According to current measurements Q rn + J1a ~ 1, 0.25 < tt m < 
0.3. The spectrum in i?(z)-coordinate for the SDSS catalogue after subtraction the background using the 
procedure described above, and applying to the residual the Rayleigh criterion Q is shown in Fig. [8] Again, 
near the values z* = 0.06 and 0.1 there are the maxima in the spectrum, but with a low confidence level. 
Moreover, it seems that the maximum at = 0.06 is undistinguishable from the near maxima, and the 
peak at z* w 0.1 is only slightly higher than its neighbours. In general the plots in Fig. [7] and Fig. [5] do not 
differ enough from each other. 

The point here is that the functions R(z) and ln(l + z) closely approximate each other. Both functions 
are shown in Fig. [9] and in the interval for < z < 7 they differ not more than 10%, i.e. exactly in the 
interval in which drop the quasars in the SDSS catalogue. For large z the functions behave in different ways: 
when z — > oo R(z) tends to a horizontal asymptote (for O m = 0.28 and ft\ — 0.72 this value is 3.3988), 
while ln(l + z) tends to infinity. 

Thus, in this case we also cannot confirm that in the quasar distribution over redshift exists a periodic 
component. 
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Figure 6: Quasar distribution over redshift in the SDSS catalogue after subrtacion the background compo- 
nent. 
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Figure 7: Quasar distribution spectrum K\{T) in the SDSS catalogue after subtaction the background 
component. 
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Figure 8: Quasar distribution spectrum Ki(T) according to the SDSS catalogue. The absciss axis corre- 
sponds to the period, espressed in the units of geodesic cosmological distance R{z) from (jHJ. 




Figure 9: The plots of functions R(z) from © for fl m = 0.28, fl A = 0.72 and ln(l + z). When z -» oo the 
function R(z) tends to a horizontal asymptote, but ln(l + z) goes to infinity. The lines cross at z = 6.60517. 
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4 Discission and conclusions 



From the analysis above one can draw a conclusion that the reliable extraction of a regular periodic com- 
ponent in the the quasar distribution over redshift is failed. However, the rippling of the quasar density, 
probably, esceeds the statitical errors. 

Most likely we deal with a cellular structure. The individual cell walls, appearing on the line of sight make 
their contribution to the numbre of quasars at fixed z. The quasar z-distribution structure of such a kind 
is not already stochastic, though it retains many properties of the random distribution. In particular, the 
dispersion may appear greater than theoretical, because an average amount of quasars inside the "bubble" 
and in its walls differs significantly from the averaged numbre of the quasars in the unit volume. 

Except that at present time the SDSS catalogue covers less than a quarter of the celestial sphere, i.e. 
the quasar distribution over the celestial sphere in it is essentially non-homogeneous. Further development 
of the catalogue should smooth all irregularities in z-distribution. As a result the spectral features in the 
distribution will be revealed at gradually decreasing confidence level. 

Summarizing our discussion, one can draw a conclusion that the periodic component in quasar z- 
distribution is absent at high confidence level. 
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